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Abstract 

t-H , 

C^l The parity- violating asymmetry arising from inelastic electron- nucleon scattering at backward angle (~ 95°) near the 
' # J °A(1232) resonance has been measured using a hydrogen target. From this asymmetry, we extracted the axial transition 
-form factor G^ A , a function of the axial Adler form factors Cf. Though G^ A has been previously studied using charged 
[current reactions, this is the first measurement of the weak neutral current excitation of the A using a proton target. 
■ -For Q 2 = 0.34 (GeV/c) 2 and W = 1.18 GeV, the asymmetry was measured to be -33.4 ± (5.3) stat ± (5.1)^ ppm. The 
Value of G^ A determined from the hydrogen asymmetry was — 0.05± (0.35) s t a t± (0.34) S1)S ± (0.06) theory. These findings 
agree within errors with theoretical predictions for both the total asymmetry and the form factor. In addition to the 
hydrogen measurement, the asymmetry was measured at the same kinematics using a deuterium target. The asymmetry 
for deuterium was determined to be —43.6 ± (lA.6) sta t ± (6.2) sys ppm. 

Keywords: electroproduction, Delta, axial, form factor, inelastic, asymmetry 
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1. Introduction 

A measurement of the parity-violating (PV) asymme- 
try in inelastic electron-nucleon scattering near the A res- 
onance has been performed as part of the G° experiment 
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P|. The study of the PV asymmetry in A electroproduc- 
tion using a neutral weak probe was first considered by 
Cahn and Gilman [2] as a test of the Standard Model, due 
to the strong dependence of the predicted asymmetry on 
the weak mixing angle, Ow- However, it is now clear that 
uncertainties in the theoretical predictions of the asym- 
metry arising from hadron structure and radiative correc- 
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tions prevent a precision Standard Model test using this 
reaction. Instead, the inelastic asymmetry can be used 
to study other physics topics, including the axial N — » A 
transition. In this paper we report the first measurement 
of this neutral current excitation of the A with hydrogen 
and deuterium targets. The axial transition form factor, 
G^ A , a linear combination of the axial Adler form factors, 
was determined from the results on hydrogen. 

A production occurs in both the charged current (CC) 
and neutral current (NC) channels of the weak interaction. 
In the constituent-quark model for the nucleon, two of the 
quark spins are aligned while the third is anti-aligned. The 
A results from interactions that flip the spin of the anti- 
aligned quark, leading to a total A spin of J = |. In 
CC reactions, there is a quark flavor change in addition 
to the spin flip. Theoretical predictions for the PV asym- 
metry rely on isospin symmetry and the assumption that 
the CC form factors will not differ significantly from their 
NC analogues. A precise measurement of the inelastic PV 
asymmetry could confirm that these assumptions are valid. 

The axial vector response in the N — > A transition has 
been studied theoretically using chiral constituent quark 
models Q, chiralperturbation theory 0, 0|, light cone 
QCD sum rules [(? . and lattice gauge theory 0, H, 0]. 
Available data are sparse, however, and the only data from 
nucleon targets come from CC neutrino experiments using 
bubble chambers [l(J 11. 12|. This subject is also topical 
because lack of knowledge of G^ A is one of the dominant 
uncertainties in understanding neutral current pion pro- 
duction [l3j |. an important background in v e appearance 
neutrino oscillation experiments. 

2. Theoretical Background 

In the scattering of longitudinally polarized electrons, 
the interference between 7 and Z° exchange amplitudes 
leads to a parity- violating dependence of the cross section 
on the electron helicity. One can form an experimental 
asymmetry 

^ _ cLor - da L 
dan + daL 

where dan and da^ are the differential cross sections for 
scattering of right- and left-handed electrons, respectively. 
A measurement of this asymmetry then allows study of 
the neutral-current contribution to the reaction. 

The PV asymmetry in electron-proton scattering near 
the A resonance can be written 14 1 



A- , - -A 
= j4 x + 



A 71 ' 
A (l) 

A 2 + A 3 , 



A (2) 



A (3) 



(2) 



where the use of the ir superscript indicates single pion 
production and ^4° is defined as 



A° = 



G F Q 2 



(3) 



where Gf is the Fermi constant, a is the fine structure 
constant, and Q 2 is the negative four- momentum transfer 
squared (Q 2 = —q 2 )- A evaluates to —61 ppm at the 
experimental kinematics (see Section [3]). The three A^ 
terms represent a decomposition of the hadronic response 
into resonant vector, non-resonant vector, and axial vec- 
tor contributions, respectively. Here, "vector" and "axial 
vector" refer to the Lorentz structure of the hadronic cur- 
rent in the scattering. In the discussion that follows, we 
present our formalism for Ai ne i which combines formalisms 
available in the literature to develop our method for the 
determination of the axial response. The present formal- 
ism was primarily derived from the work in Refs 14 1, (l5j 
and [161 ]. 

The resonant vector hadron term, A 7 ^, is the domi- 
nant term in the asymmetry and the only one that is not 
dependent on hadronic structure. It depends only on the 
weak mixing angle, Q\y. A?V contains the full contribu- 
tion of the resonant vector current at the hadronic vertex 
to the asymmetry and can be written 



14] 



A (i) 



2(1 - 2 sin 2 6 W ) 



(4) 



where g\ is the axial vector coupling of the electron to 
the Z boson, which is equal to 1 in the Standard Model, 
and £y- =1 is the isovector hadron coupling to the vector 
Z, which is 2(1 — 2 sin 2 Ow) in the Standard Model. Us- 
ing the world value sin 2 W = 0.2353 ± 0.00013 Ii7). ATj, 
is computed to be 1.06, leading to an A\ = —32.2 ppm 
contribution to the asymmetry. 

Two different formalisms were applied to the terms 
AJ 2 j and A^ . The first is a phenomenological approach, 
outlined by Musolf et al. [l4| • The second uses the dynam- 
ical model of Matsui, Sato, and Lee [l8[, and is described 
later. 

In the phenomenological approach, the non-resonant 
vector term, A 7 ^, is written as a sum of longitudinal, 
transverse magnetic, and transverse electric multipoles us- 
ing an isospin decomposition 14J. We computed the sum 
through I = 2 using multipoles from MAID2007 0. Fur- 
ther details on our implementation are available elsewhere 
(20l |. The theoretical value of A^ was determined to be 
0.018, leading to an asymmetry of A2 = —0.55 ppm. An 
uncertainty of o^' 1 = 0.72 ppm is applied to account for 
estimates made in the model. 

The axial term, A 7 ^, can be written in terms of F(Q 2 , s), 
which contains both axial and electromagnetic form fac- 
tors, 



A 1 
A (3) 



2(1- Asm 2 9 w )F{Q 2 ,s) , 



(5) 



' _1 have 



where s is the Mandelstam invariant and gy and £ A 
been replaced with their respective Standard Model tree 
level values of g v = (—1 + 4 sin 2 9 W ) and £l =1 = ~ 2. In 



Ref. [14:], A 7 ^ is defined as the total axial contribution, 



2 



including both resonant and non-resonant terms. Here, 
we neglect any non-resonant axial contributions (hence the 
use of "sa" in Equation [5]), because theoretical studies 21 1 
22J indicate that these are small. 

The function F(Q 2 ,s) can be written as a product of 
two functions of form factors, 



F(Q 2 ,s) 



E + E' 



M 



(0) 



where H EM (Q 2 ,9) and G A rA (Q 2 ) are linear combinations 
of electromagnetic (C' ; 7 ) and axial (C A ) form factors, re- 
spectively, and M is the nucleon mass. The explicit func- 
tional form of Gfj A (Q 2 ) will be presented below, while 
H EM (Q 2 ,0) is described in Ref. [20|. 

In order to calculate a theoretical asymmetry, it is nec- 
essary to compute the form factors and C A . One con- 
venient way to express the Q 2 dependence of the form 
factors is through the use of dipole forms. In this nota- 
tion, referred to as the Adler parameterization 23] 24], the 
form factors are expressed as 



CUQ 2 ) = C?{Q)G V d(Q 2 ) , 
C A (Q 2 )^C A (0)G A (Q 2 )C A (Q 2 ) 



(7) 
(8) 



where the functions G^j A (Q 2 ) are dipole form factors de- 
fined as 



G V » A (Q 2 ) = 



Q 2 



(9) 



The parameters My, a are the vector (V) and axial (A) 
dipole masses, which have been determined from fits to 
existing data. The current world values for these masses 
are M v = 0.84 GeV and M A = 1.03 ± 0.02 GeV 0. 

The function £ A is used to give additional structure to 
the Q 2 dependence of the axial response and is written 



H A (Q 2 



a'Q 2 



Q 2 



(10) 



with the parameters a 1 and b' determined from a fit to CC 
neutrino data performed by Schreiner and von Hippel 27 1. 
For the Adler model, a 1 was found to be —1.2 and b' was 2 
(GeV/c) 2 . These results hold only for Q 2 < 0.5 (GeV/c) 2 , 
which covers the present experimental kinematics. 

The values for d (0) are determined from fits to charged 
current data and are fit-dependent. In this work, the Adler 
values of these coefficients, as quoted by Nath [15j, were 
used. They are 



C 3 » = 0, 
Cf (0) = -0.35 
C#(0) = 1.20 . 



C 1 ^) = 1-85 , 
C2(0) = -0.89 



(11) 



Note that not all the Ci 's contribute to the asymmetry; Cg 
vanishes due to CVC and C A vanishes if we neglect the 



electron mass. From the coefficients above, we can also 
neglect C A . Further, the photo- and electroproduction 
data can be fit usin g th e assumption that Cj = and that 



M si7 
' M+M A 3 



28j . Thus, only the i — 3,4 terms of 



CI = 

the electromagnetic and the i = 4,5 terms of the axial form 
factors contribute to the asymmetry. The axial transition 
form factor expressed in terms of the Adler form factors is 
then 



Cna(Q 2 



\[M 2 



Ml + Q 2 ]C£{Q 2 ) - M 2 C A {Q 2 



(12) 



where M A is the mass of the A(1232). 

Using this formalism leads to the value A3 = —1.8 ppm 
for the resonant axial-vector hadron contribution to the 
asymmetry. We assigned an uncertainty of 03 = 0.65 ppm 
to this value to account for the neglect of non-resonant 
contributions and uncertainties related to the model and 
fits used to compute H EAI . Summing the three individ- 
ual components leads to a total theoretical asymmetry of 
A th = —34.6 ±1.0 ppm at the present kinematics. 

Electroweak radiative corrections (EWRC) have been 
applied to the theoretical asymmetry. We distinguish be- 
tween 'one-quark' and 'multi-quark' corrections, adopting 
the terminology of Zhu et al. 16]. One-quark radiative 
corrections are those in which the electron interacts only 
with a single quark in the nucleon, and are calculable with 
sufficient precision. In contrast, multi-quark corrections 
for the axial response are largely model dependent. 

One-quark EWRC were applied to each A^ term in- 
dependently using the corrections reported by the Parti- 
cle Data Group [17j for the MS renormalization scheme. 
These corrections have been included in the theoretical 
asymmetries quoted previously. The one-quark EWRC are 
small (< 1.5%) for both of the vector hadron terms. For 
the axial term, however, these EWRC are large, leading 
to a 58% reduction in A3. 
Zhu 



et 



161 ] have studied the multi-quark axial 



EWRC and their possible impact on the PV asymmetries. 
Specifically, they highlighted two effects present at the PV 
7./VA vertex that may contribute to the present measure- 
ment. The first effect is an anapole moment analogous to 
that previously seen in elastic scattering 29]. The second, 
referred to as the Siegert term, is of interest because it 
could lead to a non-zero asymmetry at Q 2 = 0. In a sepa- 
rate analysis, as part of the G° experiment, we have taken 
pion photoproduction data on a deuterium target at very 
low Q 2 [30j. These data allow us to bound the impact of 
the Siegert term at the present kinematics to |^4si eS ert| < 
0.15 ppm. Given the large overall theoretical uncertainty 
in these effects, multi-quark EWRC are not applied in the 
present analysis. 

The second asymmetry formalism is that of Matsui, 
Sato, and Lee who have developed a dynamical model 
of pion electroproduction near the A resonance (l8| and 
performed a calculation of the inelastic asymmetry at the 
present kinematics. As in the phenomenological approach, 



3 



Matsui et al. write the asymmetry as a sum of reso- 
nant vector, non-resonant vector, and axial vector hadron 
pieces. In their notation, Ai ne i is given by 



(2 -4 sin 2 W ) + A\ 



(13) 



where (2 — 4 sin 9w) is identical to as defined in 

Equation [4]- The two remaining terms, Ay and A_a, are 
equivalent to AJ 2 ^ and A^, respectively. However, the 
formalism used to calculate these terms differs from that 
which was presented previously. The differences in these 
two formalisms will be discussed here; more detailed in- 
formation on the dynamical model is available in Refs 18] 
and [H|. 

For the term Ay, Matsui et al. derive an expression 
in terms of structure functions analogous to that used to 
determine the resonant form of A^ as presented in Ref 



15| . This allows them to determine the contribution of the 
non-resonant vector hadron reactions to the asymmetry 
using their dynamical model rather than through the use 
of phenomenological multipoles. 

For the axial term, the definitions of Aa and A^ in 
terms of structure functions are the same. Where Matsui 
et al. differ is in their parameterization of the form fac- 
tors. As before, a dipole form is used for both the vector 
and the axial vector form factors. However, the additional 
Q 2 parameterization present in the function £ takes on 
an exponential form rather than that of Equation [TUJ Ex- 
plicitly, £ is given by 



C A (Q 2 ) = (l + aQ 2 )e 



-bQ 2 



(14) 



where a = 0.154 (GeV/cT 2 and b = 0.166 (GcV/c)- 2 
were determined by fits to CC neutrino data. The re- 
sulting theoretical values for and A3 computed with 
the phenomenological approach and the dynamical model 
each differ by < 1 ppm. 

3. Experiment 

The G° experiment was performed using a beam of 
longitudinally polarized electrons provided by the acceler- 
ator at the Thomas Jefferson National Accelerator Facility. 
The electrons were scattered from a 20 cm long unpolar- 
ized liquid hydrogen or deuterium target. The G° exper- 
iment ran in two phases, collecting data at both forward 
and backward scattering angles. The inelastic measure- 
ment was performed in the backward angle configuration. 
Since the primary concern of the G° collaboration was the 
elastic electron scattering measurement, the experimen- 
tal apparatus was optimized for elastic kinematics. For 
the backward angle measurement, the G° spectrometer 
consisted of a toroidal magnet, scintillation detectors and 
aerogel Cherenkov counters. Measurements were taken at 
beam energies of 687 MeV and 362 MeV, though only the 
higher of these energies is relevent to the inelastic mea- 
surement. Detailed descriptions of each of the components 



of the G° experimental apparatus, and their use in both 
configurations, are given elsewhere (32| . 

Electrons that scattered from the target were bent in 
a magnetic field and passed through a collimator system 
before entering the detector system. The collimators de- 
fined the experimental acceptance, leading to an effective 
scattering angle of ~ 95° for inelastic events and ~ 108° 
for elastic events. The detector system was segmented into 
octants arranged symmetrically around the beamline with 
each detector octant covering one of the eight gaps between 
adjacent magnet coils. Each detector octant consisted of 
two sets of plastic scintillators and an aerogel Cherenkov 
detector. The scintillators, labeled Focal Plane Detec- 
tors (FPDs) and Cryostat Exit Detectors (CEDs), were 
used for a rough tracking of the scattered electron's path. 
This led to a two dimensional CED-FPD detector space 
which allowed for a kinematic separation between elasti- 
cally and inelastically scattered electrons. The Cherenkov 
detector (CER) was used to distinguish between electrons 
and pions. A DAQ system consisting of specially-designed 
logic boards counted coincidences of these detectors, with 
CED-FPD-CER coincidences counted as electron events 
and CED-FPD-CER counted as pion events. Thus, elec- 
trons and pions were counted side-by-side and their rates 
recorded separately by scalers. 

The helicity of the beam was flipped at a rate of 30 
Hz, resulting in a series of 1/30 s segments of common 
helicity called macropulses (MPSs). The helicity pattern 
was generated as a collection of four MPSs, referred to as a 
quartet. The use of quartets, coupled with the fast helicity 
reversal, cancels linear drifts that can affect the asymme- 
try. The sequence of the helicity reversal for each quartet 

was chosen to be either H h or — h H — depending on 

a randomly generated initial MPS. The coincidence count 
for each CED-FPD pair was recorded for a single MPS, 
then normalized to the beam current to create an MPS 
yield. The asymmetry was then computed for each quar- 
tet by combining the detector yields for the positive and 
negative helicity MPSs according to 



*-qrt 



(15) 



Figure [T] shows the octant-averaged electron yields for the 
hydrogen data for each CED-FPD coincidence cell. Co- 
incidence cells with common kinematic acceptance were 
grouped into loci. The inelastic and elastic loci are indi- 
cated on the figure. The measured asymmetry for a given 
locus is taken as the average of the asymmetries in each 
locus cell, weighted by statistics. 

The average kinematics for the inelastic measurement 
were determined through the use of a GEANT3 simulation 
[33^ . The detector acceptance is defined by the collima- 
tors that are part of the spectrometer. For events in the 
inelastic locus, the range of accepted scattering angles was 
85° < 9 < 105° and the average 9 = 95°. This results in 
a Q 2 range of 0.25 < Q 2 < 0.5 (GeV/c) 2 with an average 
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O 9 




3 4 5 6 7 



9 10 11 12 13 14 15 16 



Figure 1: Octant averaged electron yields for the high energy hy- 
drogen run period. On the j/-axis is CED number and the x-axis is 
FPD number, and each block represents a CED-FPD coincidence. 
The color scale indicates cell yield in Hz//xA. The collection of cells 
outlined in with a solid (dashed) line represents the inelastic (elastic) 
locus. Note that since FPDs 1 and 2 were not used in the backward 
angle configuration, the FPD numbering begins with 3. 



of Q 2 = 0.34 (GeV/c) 2 . The range of invariant mass, W, 
covers the entire region from the pion threshold to just 
past the peak of the A, with a range of 1.07 < W < 1.26 
GeV. Most of the events were on the low W side of the 
peak, leading to an average just below the resonance at W 
= 1.18 GeV. 



4. Data Analysis 

In order to determine the physics asymmetry, several 
corrections were applied to the raw asymmetry An overview 
of the corrections will be provided in this section, while 
more detailed descriptions of the data analysis related to 
beam and instrumentation [34[ > and the inelastic result [20| 
are available elsewhere. A summary of the corrections ap- 
plied and their contributions to the systematic uncertainty 
is given in Table O 

There were three rate-dependent effects that needed to 
be considered. Corrections were applied to the measured 
yield to account for two of these: dead time in the detec- 
tors and electronics, and randomly triggered Cherenkov 
events. Typical dead times were ^30 ns. Applying these 
corrections resulted in an increase in the magnitude of the 
asymmetry for both targets. For hydrogen, there was a 2 
ppm, or 10%, increase, while for the deuterium the change 
was much larger at 12.6 ppm, or nearly 90% of the un- 
corrected value. This large shift in the asymmetry for 
the deuterium data was due to the Cherenkov corrections, 
which were more significant because of high rates of 7r~ . 
Error analysis showed that for each data set the two cor- 
rections applied provided a negligible contribution to the 
uncertainty. 

The third rate-dependent effect, random CED-FPD co- 
incidences, was not included in the corrections applied. 



This effect was instead treated as an uncertainty. For each 
data set, an upper bound on the residual false asymmetry 
due to the presence of CED-FPD randoms was estimated 
using information obtained from the elastic electron lo- 
cus, where rate-correlated effects were more pronounced. 
The bounds determined from this analysis, 0.16 ppm for 
hydrogen and 1.2 ppm for deuterium, were assigned as un- 
certainties for the rate corrections. 

Corrections were also applied to account for any false 
asymmetry arising from helicity-correlated beam proper- 
ties, including angle, position and energy of the electron 
beam. Because of the high quality of beam that was pro- 
vided by the Jefferson Lab accelerator, helicity-correlated 
effects were negligible. As a result, the impact of these 
effects was also small (Af a i se < 0.3 ppm). A conservative 
100% uncertainty was assigned for this correction. 

The beam polarization was determined using a M0ller 
polarimeter located upstream of the target [35J. Mea- 
surements taken periodically throughout the experimen- 
tal run found the polarization to be steady at P = 85.8 
% ± (0.07) sta4 ± (\A) sys . Though the beam polarization 
was nominally longitudinal, there was a small transverse 
component present. However, the symmetrical nature of 
the G° spectrometer meant this component's impact on 
the asymmetry was negligible and no correction was nec- 
essary. Instead, using information determined through a 
measurement of the asymmetry with transverse beam and 
an estimate of detector misalignment, an upper bound on 
the false asymmetry was determined for each target and 
treated as an uncertainty (See Table H]) 

The most significant correction applied to the inelas- 
tic asymmetry was the background correction. Since the 
spectrometer was optimized for elastic, not inelastic, scat- 
tering, the yield in the inelastic region of the matrix con- 
tained a high percentage of background. In order to cor- 
rect for backgrounds, a procedure was developed 20] that 
made use of both background measurements and simulated 
yields to determine the fractional contribution, or dilution 
factor (jf 9 ), of each process. The dilution factors were 
then used to subtract the background asymmetries (A^ 9 ) 
from the cell average asymmetry ( A meas ) according to 



A 



in el 



A — V f bg A bg 

-n-meas J j -^j 



i-E/, 



bg 



(16) 



In any given CED-FPD cell, there were up to five major 
processes contributing to the total yield and average asym- 
metry: electrons scattered elastically from the target liq- 
uid, electrons scattered inelastically from the target liquid, 
electrons scattered elastically or inelastically from the alu- 
minum windows of the target cell, ir° decay and misiden- 
tified 7r~. Note that the ir~ contribution here differs from 
that which is removed by the rate corrections, with the 
contamination resulting both from rate-independent ef- 
fects causing false Cherenkov triggers and from <5-ray pro- 
duction. The 7T° contribution to the background is due to 
electrons that are emitted through secondary processes fol- 
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lowing pion decays, primarily electron-positron pairs from 
decay photons interacting in the shielding. 

The yields due to scattering from the target windows 
and from pion contamination were determined using data 
from special measurements made during the experimental 
run. For the three remaining processes, simulation was 
used to model the individual yield distributions across the 
detector acceptance while a fitting routine determined the 
appropriate normalizations. Before applying the fit, the 
target window and n~ yields were subtracted from the 
total yield leading to a reduced yield, Yr, consisting of 
elastic, inelastic and tt° decay yields. Yr was then fit as a 
function of FPD for each CED, or across each row in the 
yield matrix shown in Fig. [TJ by assigning scale factors to 
the simulated yields and allowing them to vary indepen- 
dently. Since our symmetrical spectrometer implies con- 
stant yields across the detector octants, all octants were 
fit simultaneously to find a single scale factor that rep- 
resented an octant average for the CED. Once the scale 
factors were determined, a second fit of scale factor as a 
function of CED was performed for each process to ensure 
that the scale factors varied smoothly across the matrix. 
The yield contribution for a given process in a given cell 
was then defined as the simulated cell yield multiplied by 
the fitted scale factor for that cell. Details on the meth- 
ods used to determine the cell-by-cell yield contribution 
from each process are given in [20[ • Figure [5] shows the 
full background separation for a typical CED in a typical 
octant for the hydrogen data. 




FPD 

Figure 2: Contributions of the different processes present in CED 
4, Octant 3, for the hydrogen data. The data are given as points 
and the solid curve represents the total yield as determined through 
background analysis. The individual processes are shown as curves: 
inelastic scattering (dot), elastic scattering (dash), n° decay (dash- 
dot), scattering from Al target windows (long dash). The vertical 
line is included to highlight cells located in the elastic and inelas- 
tic loci, as indicated by the labels. Note that the yields determined 
through background analysis are drawn as curves rather than indi- 
vidual points to make the plot more clear. The shape of the curves 
comes only from connecting the points together and has no other 
significance. 

Measured asymmetries determined from 67° data were 
used for each of the background processes. An asymme- 
try of ± 3 ppm, based on our measured ir~ asymmetry 



and experimental uncertainty, was used for both the ir° 
decay and tt~ backgrounds. Our measured elastic elec- 
tron asymmetry (36j was used with a scaling, determined 
through simulation, to account for the fact that the back- 
ground in the inelastic locus is mainly from the elastic 
radiative tail, which has a lower average Q 2 than our elas- 
tic data. The scaling was computed cell-by-cell, with those 
closest to the elastic locus having the smallest scale factors. 
The locus average scaling was about 30% for hydrogen and 
over 50% for deuterium. For the empty target background, 
the dominant process was A-electroproduction from alu- 
minum. Although we were able to use data to determine 
the yields, the data did not have sufficient precision for an 
asymmetry measurement. Instead, the aluminum asym- 
metry was approximated using the measured deuterium 
asymmetry. This asymmetry differs slightly from the true 
aluminum asymmetry due to the differing proton to neu- 
tron ratio between the two nuclei, however, this difference 
is small compared to the error in the measurement. Note 
that it was not necessary to remove the empty target back- 
ground from the deuterium data since the asymmetries 
were considered equal. 

Tablc[T]summarizcs the dilution factors and their corre- 
sponding background asymmetries. The total background 
fraction in the inelastic locus was (52.5 ±4)% for the hy- 
drogen data and (64.9 ± 5)% for the deuterium data. For 
both targets, the radiative tail of the elastic electrons was 
the dominant background, contributing 25% of the hydro- 
gen yield and 30% of the deuterium yield. The higher 
total background fraction in deuterium was due both to 
the presence of 7r~ contamination and to the widening of 
the elastic peak due to Fermi motion in the deuterium nu- 
cleus. Subtraction of the backgrounds results in a —7.3 
ppm change in the hydrogen asymmetry and a —12.6 ppm 
change in the deuterium asymmetry. The background cor- 
rection also results in an approximately 5 ppm systematic 
error due to the uncertainties from both the measured 
asymmetries and the dilution factors. 



Hydrogen 687 MeV 



Process 


fbg 




Abg 


Elastic 


25.7 ± 0.4 % 


-14.5 ± 0.89 


ppm 


Empty Target 


15.6 ± 0.8 % 


-43.6 ± 16 


ppm 


7T° Decay 


13.6 ± 3.6 % 


± 3 


ppm 




negligible 


± 3 


ppm 


Deuterium 687 MeV 


Process 


fbg 




Abg 


Elastic 


31.0 ± 0.3 % 


-31.4 ± 2.2 


ppm 


Empty Target 


9.0 ± 0.5 % 






7T° Decay 


11.3 ± 3.2 % 


± 3 


ppm 


7T~ 


11.1 ± 3.3% 


± 3 


ppm 



Table 1: Summary of the dilution factors and background asymme- 
tries. 



For the hydrogen asymmetry, two additional correc- 
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tions were needed before comparing our results to theory. 
The first of these accounts for electromagnetic radiative 
effects by using simulation to compute corrections accord- 
ing to the procedure of Mo and Tsai 37 1 [38| . The to- 



Hydrogen 687 MeV 



dA (ppm) a c s ° r s (ppm) 



tal correction was determined by comparing the simulated 
locus-average value of the asymmetry with and without 
radiative corrections included and was found to be (1.17 
± 0.6)%. The uncertainty arises from the model chosen for 
the inelastic cross section. We then accounted for bias due 
to acceptance by comparing the locus average asymmetry, 
(A(Q 2 , W)), to the asymmetry at the average kinematics, 
A{(Q 2 ), (W)), and found this to be a (-1.6 ± 0.6)% effect. 
The uncertainty here was determined by comparing the 
asymmetries computed with the two formalisms discussed 
in Section [2] 

Table [2] summarizes the corrections applied, the sys- 
tematic uncertainty assigned for each, and their impact 
on the asymmetry. Backgrounds constituted the largest 
correction to the asymmetry and the largest contribution 
to the uncertainty. Note that, in addition to its impact 
on systematic uncertainty, the background correction led 
to an increase in the statistical uncertainty as background 
events were removed. Additionally, since the hydrogen re- 
sult made use of the deuterium result for the aluminum 
asymmetry, the systematic uncertainty and the large sta- 
tistical error present in the deuterium asymmetry were fac- 
tored in to the systematic uncertainty for hydrogen. The 
second largest impact was that of the rate corrections, with 
the larger contribution to the deuterium asymmetry due 
to generally higher scattering rates and greater contami- 
nation caused by high pion rates. The total uncertainty, 
including both systematic and statistical errors, is 7.4 ppm 
for the hydrogen asymmetry and 16 ppm for deuterium. 

5. Results and Conclusion 

Upon applying the corrections outlined in the previous 
section, we arrive at the final measured asymmetries for 
the hydrogen and deuterium data. These asymmetries are 



Ai nel =-33.4 ±(5.3) stat ± (5.1) 
Ag el = -43.6 ±(U.6) stat ± (6.2) 



sys- 



(17) 
(18) 



We found the deuterium asymmetry to be consistent with 
the hydrogen asymmetry within errors as would be ex- 
pected. As discussed previously, due to the lack of a model 
for the deuterium asymmetry, no further information was 
extracted from Af nel . For the hydrogen asymmetry, how- 
ever, theoretical input can be used to determine the form 
factor of interest in this measurement. 

The measured G^ A can be extracted from the asym- 
metry by first using the theoretical values of Ai and A2 
discussed in Section [5] to isolate the axial response, A3. 
Subtracting off the vector portions of A ine i leads to A3 = 
-0.69 ± (5.3) stat ± (5.1) sys ± (0.7) t ft ppm, where the sub- 
script "th" denotes the theoretical uncertainty. The theo- 
retical value for the resonant contribution calculated using 



Rate Corrections 


-2.17 


0.16 


Helicity Cor. Beam 


-0.16 


0.16 


Beam Polarization 


-3.91 


0.36 


Transverse Polarization 




0.03 


Backgrounds 


-7.33 


4.93 


EM Radiative Effects 


-0.39 


0.20 


Acceptance Averaging 


-0.55 


0.20 


Total a sys 




5.11 


Total OStat 




5.30 


Deuterium 


687 MeV 






dA (ppm) 


VsZ (PP m ) 


Rate Corrections 


-12.6 


1.20 


Helicity Cor. Beam 


+0.25 


0.25 


Beam Polarization 


-4.66 


0.43 


Transverse Polarization 




0.02 


Backgrounds 


-12.5 


5.52 


Total a sys 




6.23 


Total a s tat 




14.6 



Table 2: Systematic uncertainty (0^°^) due to the corrections applied 
and each correction's impact (dA) on the average asymmetry. The 
total statistical uncertainty, <7 s tat is given for comparison purposes. 



the formalism of Musolf et al. is A3 = —1.8 ppm, while 
the Matsui et al. formalism leads to A3 = —1.7 ppm. The 
measured value for A3 is consistent within errors with both 
of these values. While these results indicate that there are 
no major deficiencies in the theory, or in the assumption 
that the non-resonant contribution can be neglected, there 
is insufficient precision in the present measurement to se- 
lect one prediction over the other. 

From this value of A3, G^ A was found to be —0.05 ± 
(0.35) stQt ± (0.M) sys ± (0.06) th - Here the full theoretical 
uncertainty of 1.0 ppm has been included. The resulting 



G^- A is consistent, within errors, with the theoretical value 
of -0.196. 

Figure [3] shows a comparison of the measured hydrogen 
asymmetry to the theoretical values computed using the 
definitions in Section [5J The statistical, systematic and 
theoretical uncertainties have been summed in quadrature 
to yield a single total uncertainty. The measured asym- 
metry, Ai ne i is consistent with the total theoretical asym- 
metry, A t h — A\ + A2 + A3, within errors. 

In summary, we have measured the parity-violating 
asymmetry in inelastic electron-proton scattering and found 
excellent agreement with theoretical expectations based on 
the formalisms of Musolf et al. [l4| and Matsui et al. [l8| • 
In addition, the axial transition form factor, G^ A , has 
been determined to be consistent with the theoretical for- 
malism presented by Musolf et al. 14) using the Adlcr 



parameterization. This represents the first neutral current 
measurement of the axial N— > A response from a nucleon 
target. 
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Figure 3: Measured hydrogen asymmetry (point) plotted with the 
total theoretical asymmetry, A t /j, computed using the Musolf et al. 
formalism (line). The asymmetry components Ai (long dash), A 2 
(dot) and A3 (dash-dot) are also shown. Inset: Comparison of the 
measured asymmetry to the total A t / X computed using the formalisms 
of Musolf et al. (line) and Matsui et al. (dash). 
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